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Abstract — In this paper, coded slotted ALOHA (CSA) is 
introduced as a powerful random access scheme to the MAC 
frame. In CSA, the burst a generic user wishes to transmit in 
the MAC frame is first split into segments, and these segments 
are then encoded through a local a packet-oriented code prior to 
transmission. On the receiver side, iterative interference cancel- 
lation combined with decoding of the local code is performed 
to recover from collisions. The new scheme generalizes the 
previously proposed irregular repetition slotted ALOHA (IRSA) 
technique, based on a simple repetition of the users' bursts. An 
interpretation of the CSA interference cancellation process as an 
iterative erasure decoding process over a sparse bipartite graph 
is identified, and the corresponding density evolution equations 
derived. Based on these equations, asymptotically optimal CSA 
schemes are designed for several rates and their performance 
for a finite number of users investigated through simulation and 
compared to IRSA competitors. Throughputs as high as 0.8 are 
demonstrated. The new scheme turns out to be a good candidate 
in contexts where power efficiency is required. 



I. Introduction 

Although demand assignment multiple access (DAMA) 
medium access control (MAC) protocols guarantee an efficient 
usage of the available bandwidth MAC random access 
schemes remain an appealing and popular solution for wireless 
networks. Among them, slotted ALOHA (SA) Q-Jl] is 
currently adopted as the initial access scheme in both cellular 
terrestrial and satellite communication networks |@]. In llS] an 
improvement to SA was proposed, namely, diversity slotted 
ALOHA (DSA). In DSA, each packet (also called burst) is 
transmitted twice over the MAC frame, which provides a slight 
throughput gain over SA. As a drawback, for the same peak 
transmission power of the SA scheme, the average transmitted 
power of DSA is doubled. 

A more effective use of the burst repetition is provided by 
contention resolution diversity slotted ALOHA (CRDSA) l6l, 
whose basic idea is the adoption of interference cancellation 
(IC) to resolve collisions. More specifically, with respect to 
DSA, each of the twin replicas of a burst, transmitted within 
a MAC frame, possesses a pointer to the slot position where 
the respective copy was sent. Whenever a clean burst is 
detected and successfully decoded, the pointer is extracted 
and the interference contribution caused by the burst copy on 
the corresponding slot is removed. This procedure is iterated, 
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possibly allowing to recover the whole set of bursts transmitted 
within the same MAC frame. This results in a remarkably 
improved normalized throughput S (defined as the probability 
of successful packet transmission per time slot) which may 
reach S ~ 0.55, whereas the peak throughput for pure SA is 
5 = 1/e ~ 0.37. Further improvements can be achieved by 
exploiting the capture effect ||3|, ||7|. 

In m, in irregular repetition slotted ALOHA (IRSA) was 
introduced to provide a further throughput gain over CRDSA. 
A higher normalized throughput is achieved by IRSA by 
allowing a variable and judiciously designed repetition rate 
for each burst. As for DSA, the performance improvement 
achieved by CRDSA/IRSA has a counterpart in the increment 
of the average transmitted power Since CRDSA is a specific 
instance of IRSA, in the following we will refer in general 
to IRSA. In O it is also illustrated how the iterative burst 
recovery process on the receiver side can be represented 
via a bipartite graph and how, under the assumption of an 
ideal channel estimation and of a sufficiently large signal-to- 
noise ratio (SNR), it shares several commonalities with the 
graph representation of the erasure recovery process of modem 
channel codes on sparse graphs ifTOl . ifTTI . 

In this paper, we introduce a further generaUzation of 
IRSA, named coded slotted ALOHA (CSA). The basic idea 
of CSA is to encode (instead of simply repeat) bursts using 
local codes prior to transmission in the MAC frame and to 
combine, on the receiver side, iterative IC with decoding of 
the local codes to recover from collisions. The new scheme 
turns out to be interesting especially in contexts where power 
efficiency is required. Density evolution equations for CSA 
are derived to analyze the IC process in an asymptotic setting, 
leading to the calculation of the peak asymptotic throughput. 
Numerical results are then presented to illustrate the validity 
of the proposed asymptotic analysis and its effectiveness in the 
design of CSA access schemes for a finite number of users. 

II. System Model 

Similarly to 161, |I8|, we consider a random access scheme 
where the slots are grouped in MAC frames, all with the same 
length (in slots). We further restrict to the case where each user 
attempts one burst transmission per MAC frame. 

Consider M users, each attempting the transmission of a 
burst of time duration Tsa over a MAC frame of time duration 
Tp. Neglecting guard times, the MAC frame is composed of 
iVgA = Tp/TgA slots. In SA, each user would independently 
choose one of the A^sa slots uniformly at random and would 
attempt transmission of his burst into that slot. In IRSA each 



user would generate a certain number r of replicas of his burst, 
where r may be not the same for two different users, and 
would transmit the r replicas into r slots chosen uniformly at 
random among the available A^sA slots. 

In CSA, when a user wishes to transmit a burst of time 
duration Tsa over the MAC frame, the burst is divided into 
k information sub-bursts (also called information segments), 
each of time duration Tcsa = T^A/k. The k information 
segments are then encoded by the user via a packet-oriented bi- 
nary linear block code which generates nt encoded segments, 
each of time duration Tcsa = T'^h/k. For each transmission, 
the code to be employed is drawn randomly by the user from 
a set of ric possible codes. For ft, £ {!,..., Uc] the h\h code, 
denoted by '^h, is a (n.^,, fc, d^^^^-^) code, that is it has length nh, 
dimension k, and minimum distance cf^^^. We further impose 
that has no idle bits and fulfills d^^^^^ > 2. We assume 
that, at any transmission, each user independently chooses his 
local code according to a probability mass function (p.m.f.) 
-P= [Phl^'Li which is the same for all users. Denoting again 
by Tp the MAC frame duration, the MAC frame is composed 
of Nqsa — Tp/TcsA — kNsA slots. The rih coded segments 
are then transmitted by the user over rih slots picked uniformly 
at random. Note that IRSA may be seen as a special case of 
CSA where k = 1 and each %i is a repetition code of length 
rih, and that SA is a special case of IRSA where n/, = 1 for all 
usersl^The overall rate of CSA is defined as R = k/n, where 
n := Phnfi is the expected length of the code employed 

by the generic user. Note that AP — n/k = 1/R represents 
the increment of average power with respect to pure SA. 

It is now convenient to introduce a graph representation of 
CSA, depicted in Fig. [T] Let us consider a MAC frame com- 
posed of Nqsa slots, in which M users attempt a transmission. 
The MAC frame status can be represented by a bipartite graph, 
Q = {B, S, £), consisting of a set B of M burst nodes (one 
for each burst transmitted in the MAC frame), a set S of 
NcsA sum nodes (one for each slot), and a set £ of edges. An 
edge connects a burst node (BN) 6^ G ;B to a sum node (SN) 
Sj G iS if and only if an encoded segment associated with 
the ith burst is transmitted in the jth slot. In other words, 
BNs correspond to bursts, SNs to slots, and edges to encoded 
segments. Therefore, a burst split into k information segments 
and encoded via the code %i is represented as a BN with rih 
neighbors. Correspondingly, a slot where d replicas collide 
is represented as a SN with d connections. The number of 
edges emanating from a node is the node degree. Moreover, a 
BN where is employed during the current transmission is 
referred to as a BN of type h. 

In our analysis, we rely on three assumptions. 1) Sufficiently 
high SNR. This allows to claim that, when a segment is 
received in a clean slot, it is known at the receiver 2) 
Ideal channel estimation. Under this assumption (and the 
previous one), ideal IC is possible, allowing the recovery 
of collided bursts with a probability that is essentially one. 
3) Destructive collisions. Segments that collide in a slot are 
treated as erasures. These assumptions simplify the analysis 

'We point out that CSA may be seen as a generalization also of the schemes 
proposed in 1121 . 1131 . where no IC was used. 



A^CSA sum nodes 




M burst nodes 

Fig. 1. Graph representation of CSA. Circles are the burst nodes and represent 
the M users, squai'es are the sum nodes and represent the N slots. The degree 
of a burst node is equal to the length of the locally employed code. The 
degree of a sum node is equal to the number of collided encoded segments. 
The example is for k = 2. 



without substantially affecting the obtained results, as shown 
in Ig) and ||9] for CRDSA and IRSA, respectively. 

Each coded segment associated with a BN of type h is 
equipped with information about the relevant user and with 
a pointer to the other n/, — 1 segments Q On the receiver 
side, segments which collided in some slot with those sent 
by another user are marked as lost, so that a BN is connected 
to "known" edges and to "unknown" ones. Hence, some of its 
information bursts are known, and the others unknown. At the 
generic BN (say of type h), erasure decoding of the code "^^h 
may allow to recover some of the unknown encoded and infor- 
mation segments. It is now possible to subtract the interference 
contribution of the newly recovered encoded segments from 
the signal received in the corresponding slot. If c? — 1 segments 
that collided in a SN of degree d have been recovered by the 
corresponding BNs, the remaining segment becomes known. 
The IC process combined with local decoding at the BNs 
proceeds iteratively, i.e., cleaned segments may allow solving 
other collisions. Note that this procedure is equivalent to 
iterative decoding of a doubly-generalized low-density parity- 
check (D-GLDPC) code over the erasure channel llT4l . where 
variable nodes are generic linear block codes and check nodes 
are single parity-check (SPC) codes. 

Denoting by = A'csa the number of slots (a multiple of 



fc), the logical normalized offered traffic G is given b 



itipii 



G 



kM 



(1) 



The normalized throughput S is defined as the probability of 
successful packet transmission per time slot. For example, for 
standard SA we have S = Ge^'~'. 

Finally, we recall the definition of information function of 
a linear block code ifTSl . Let G be a generator matrix for an 
{n, k) linear block code ^. The 5th un-normalized information 
function of denoted by eg, is defined as the summation of 
the ranks over all the possible submatrices obtained selecting 
g columns (with < g < n) out of G. 

^In practical implementations, the overhead due to the inclusion of pointers 
in the segment header may be reduced by adopting more efficient techniques. 
For fixed k, one may include in the segment header the code index h together 
with a random seed, out of which it is possible to reconstruct (by a pre-defined 
pseudo-random number generator) the positions of the n;j segments. 

^^In CSA and IRSA we distinguish the logical load G from the physical load 
given by -^G = G/R and representing the average number of transmitted 
segments per slot. The logical load G provides a direct measure of the traffic 
handled by the scheme. Note that the two concepts coincide in pure SA. 



III. Density Evolution, Threshold, and Stability 

The degree distribution of the SNs from a node perspective 
is defined as 



(2) 



d>0 



where is the probabiUty that a SN has degree d. 

Let us consider a user encoding his segments through the 
code '^i'h, and allocating his Uh encoded segments into nt 
slots chosen randomly. Then, the probability that the BN 
associated with this user (say U) is connected to a SN 
A may be expressed as the ratio between the number of 
ways of connecting the rih sockets of U to the N SNs 
such that U is connected to A, to the total number of 
ways of connecting the rih sockets of U to the N SNs: 
Pr{[/ is connected to A\U uses %,} = = if- 

Therefore we have: 

Tie - 

VyW is connected to = > Ph — = — . 

h=l 

Since each user selects his slots independently of all the 
other users, the probability that a SN has degree d (that is 
the probability that the SN is chosen by d users) is given by 

*d = ( 1 ( 1 - 
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where equality (a) follows from ([T). Therefore, in the limit 
where AI (and consequently, for fixed G and k, N through 
([T]i) tends to infinity, (|2]i may be written as 



d>0 
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(4) 



Using (IDi we can now express the probability that an edge 
is connected to a SN of degree d > 1 as: 



Pd 



■^dd 



(5) 
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Therefore, the degree distribution of the SNs from an edge 
perspective is given by 

p[^) =e-i^y, ^r^, ' = exp (-^G{1 x)) (6) 



d>l ^ ' 



and p{x) = '^(x). 

For given k and G, we investigate the evolution of the 
decoding process described in Section in the asymptotic 
case where M — > cxd (and consequently N ^ oo through 

Proposition 1: Assume MAP decoding is used at each 
BN. At the ith decoding iteration, let pi^i be the average 
probability that an edge carries an erasure message^ from the 

'^This is the probability that an edge is associated with an encoded segment 
that is still unknown. 



SNs to the BNs. Consider a BN where is employed and let 
(jj-''^ be the average probability that an edge carries an erasure 
message outgoing from the BN, after MAP decoding at the 
BN. Then, we have 



nh 



t=o 



(7) 



where e^'*' is the gth unnormaUzed information function of "^h- 

The proof of Proposition [T] is omitted due to space constraints. 
Note that the proof follows exactly the same argument used in 
[16 Theorem 2] to derive the expression of the EXIT function 
of a linear block code without idle bits over the binary erasure 
channel. 

Proposition 2: Assume MAP decoding is used at each BN. 
At the ith iteration, let pi_i be the average probability that 
an edge carries an erasure message from the SNs to the BNs, 
before MAP decoding at the burst node. Let qi be the average 
probability that an edge carries an erasure message from the 
burst nodes to the SNs, after MAP decoding at the BNs. Then: 



7iW 



h=l t=0 



(8) 



Proof: For all h E {!,..., Uc}, let A/j be the probability 
that an edge is connected to a BN of type h. We have 

Phnii 



n 

The proposition is proved by observing that 



9i = X] -^^^ 



(h) 



(9) 



(10) 



h=l 



and by incorporating (ITJ and (|9} into (fTOt . ■ 

The following is a well-known result from basic density evo- 
lution on the erasure channel for irregular LDPC codes llTTI . 



Proposition 3: At the ith iteration, let q.^ be the average 
probability that an edge carries an erasure message from the 
BNs to the SNs before decoding at the SNs. Let pi be the 
average probability that an edge carries an erasure message 
from the SNs to the BNs after IC at the SNs. Then: 



Pt 



l-p{l-q^)■ 



(11) 



Incorporating (|8) into (fTTT i and recalling (|6), we obtain the 
nonlinear difference equation 



Pi-i 



\nh — l—t 



h=l t=0 

x[K-t)eltl*-(i + l)ei'l^_,]} (12) 

which expresses the evolution of the average probability 
that an edge carries an erasure message at the ith decoding 



iteration. The initial value of (fTSl i shall be set to po = 1- The 
asymptotic threshold G* of the system is defined as 



G* := sup{G >0■.p^ 



as 



oo, 



= 1} . 



The threshold G* is the supremum G such that, in the 
asymptotic setting M ^ oo, the normalized throughput S 
fulfills S = G. For all values of G < G* the offered traffic 
turns into useful throughput and therefore G* is the asymptotic 
peak throughput. 

Using standard bifurcation theory, the threshold G* is equal 
to the smallest G > such that, for some < .t < 1, (x, G) 
is a solution to the system of simultaneous equations 



f{x,G)=x 
dfjx.G) _^ 
dx 



(13) 
(14) 



where 



f{x,G) :=l-cxp{-^^P;, « 



h=l 



t=0 



'(1 - a;)""-!-*} 



(15) 



and a^") := (n, - t)^^., - {t + l)eX/_i_, ■ 



A. Stability 

Difference equations such as ( fT2] i are often used to model 
discrete dynamical systems. These systems are typically an- 
alyzed as regard to the stability of their fixed (or steady- 
state equilibrium) points. A fixed point x of xi = /(.t^_i) 
is known to be locally stable if there exists e > such 
that lim^^oo = x for all xo such that |xo — x\ < e. 
The following well-known result establishes a necessary and 
sufficient condition for local stability of a fixed point. 

Lemma 1: A fixed point x of a discrete dynamical system 
xg = f{xe-i), where / : M t— !> R is a differentiable and single- 
valued function, is locally stable if and only if \f'{x)\ < 1. 

It is readily shown that p = is a fixed point of ( fT2] i. 
corresponding to successful IC. Therefore we may apply 
Lemma [U to study its stability. We obtain the following result. 

Proposition 4 (Stability condition): For h G {l,...,7ic}, 
let "^/j be a (f^?i7 fc, d|^in) linear block code employed by 
each user with probability Ph at each transmission and aH^^ 
be the number of weight-it) codewords of "^tah- Moreover, let 
d = mmh{d\^l} and V = {h : d[^l = d}. If = 2, then the 
fixed point p = of (fTSI) is locally stable if and only if 

^ (16) 



G< 



2A2 



where A2 ~ J^hev PhA^2^ is the average number of weight-2 
codewords. Else, if d > 3, then the fixed point p = of (fT2] i 
is stable for any value of G. 

Proof: Let us define again f{x,G) as in dlSl l and let us 
denote by sjj^^ the generic (fc x 17) matrix obtained by selecting 
g columns in (any representation of) the generator matrix of 



^^h , irrespective of the order of the g columns, and by ^ , 
the summation over all (") matrices Sg'^\ We have: 



G 



(h) 
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(a) 
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h=l 
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(b) [ i^A2 if d=2 







if d > 3 



(17) 



In the previous equation list, both (a) and (b) rely on the 
hypothesis that dmin > 2 and on ifTTl Proposition 2]. ■ 

The stability condition is a necessary, but in general not 
sufficient condition for successful decoding. Note also that 
the stability condition implies 

k 



G* < 



2A2 



(18) 



that will be referred to as the stability upper bound, denoted by 
G*^. Note that in the IRSA case (k = 1) we have A2 = P2, 
where P2 is the probability to select the length-2 repetition 
code, and therefore for IRSA we obtain G* < 2^7- 

In the case where d = 2, ( fTSl l may be achieved with equality. 
Indeed, this is the case when iic = 1 and the binary linear 
block code 'rf employed by all users is a SPC code. 

Proposition 5: Let ric = 1 and the linear block code 
employed by all users be a (fc + 1, fc) SPC code. Then 

1 



G* = 



(19) 



k+1 

and ( fTsT i is achieved with equality. 

The proof is easily obtained by recasting (fT?t and (fT4t for 
the special case of SPC codes and by showing that {x, G) = 
(0, ^) is a solution to the system and that no G < ^ exists 
such that {x, G) is a solution to the system for any < x < 1. 

IV. CSA FROM Random Linear Block Codes 

So far the generic user has been assumed to encode, at each 
transmission, its k information segments via an (nh,k,d[^^^) 
binary linear block code picked randomly with p.m.f. P= 
[PjiYl'Li from an ensemble of ric candidate codes. In this 
section, we consider a slightly different situation. Specifically, 
we assume that, at each transmission, the generic user picks 
randomly a codeword length Ug > k from the ensemble 
{ni, . . . , Us^^^} with p.m.f. Q= [Qn^Ys^T '^'^'^ encodes his k 
segments through a binary (fc x n^) generator matrix generated 
uniformly at random from the set of all rank-fc (fc x n^) binary 
matrices representing [rig^k) linear block codes without idle 
bits and with minimum distance at least 2. We are interested 
in calculating the expected threshold G* for this scheme. The 



TABLE I 

Optimized PROBABILITY distributionP for IRSA schemes with rates 1/3, 2/5 and R = 1/2 and optimized probability distributionQ 

FOR CSA schemes with k = 2 AND RATES 1/3, 2/5, 1/2 AND 3/5 UNDER THE RANDOM CODE HYPOTHESIS. 



IRSA 


G* 


Gl 






(2,1) 


(3,1) 


(6,1) 








R = 


1/3 


0.554016 


0.261312 


0.184672 




0.8792 


0.9025 


R = 


2/5 


0.622412 


0.255176 


0.122412 




0.7825 


0.8033 


R = 


1/2 


1.000000 








0.5000 


0.5000 


CSAk = 2 


G* 








(3,2) 


(4,2) 


(5,2) 


(8,2) (9,2) (12,2) 






R = 


1/3 


0.088459 


0.544180 


0.121490 


0.245871 


0.8678 


0.9427 


R = 


2/5 


0.153057 


0.485086 


0.135499 


0.114235 0.112124 


0.7965 


0.8391 


R = 


1/2 




1.000000 






0.6556 


0.7500 


R = 


3/5 


0.666667 


0.333333 






0.4091 


0.4091 



advantage of a random code hypothesis is to allow to release 
the analysis from considering a specific set of ?ic codes. 

With respect to the previous case, the expression (2) of ^'(a;) 
and the expression (|6]l of p{x) remain unchanged, provided the 
definition of n is updated as n = X]s=r Qn^ns- Analogously, 
(fTTT i is not affected by the random code hypothesis. On the 
other hand, we now update ^ by replacing qi with its average 
value qi. Denoting by G(„^ j.) the ensemble of all rank-A: {k x 
Us) binary matrices representing linear block codes without 
idle bits and with minimum distance at least 2, and by Eg(„ 
the expectation operator over G(„^ j.), we have 



where Eg„_ ^ [at) 



1 ^'>"•- 



t=0 



Us-l-t 



X [{us - t)EG„.^ ,)(e„^-t) -{t + l)EG,„^_,,(g„,-i-t)] 

where again n = X]s=r Qn^ns- The expectation Eg(„^ (e^) 
may be calculated using the following result developed in 
ifTsl . where a recursive technique to calculate the functions 
J{k,n,k) and K{k,n, g,u, k) is also available. 

Proposition 6: For given positive integers n, k < n, and 

g <n, Eg,„ j^)(eg) is given by 



n\ -J^^^ K{k,n, g,u,k) 



E 



J(fc, n, k) 



(20) 



where J(k, n, k) is the number of rank-fc (k x n) binary matri- 
ces without zero columns and without independent columns, 
and where K{k,n, g,u, k) is the number of rank-fc {k x n) 
binary matrices without zero columns, without independent 
columns and such that the first g columns have rank 

The average threshold G* may be calculated by properly up- 
dating the simultaneous equations ( fTsl l and (fT4l i. Specifically, 
defining the function f{x) as 

/>, G) := 1 - cxp { - - £ g„^ J^G.. 

X -x)"--!^*} (21) 

^In this context, a column is said to be independent when deleting the 
column from the matrix does not affect the rank of the matrix. 



[{us - i)EG(„^ ,,(g„,-t) - (i + 1) 



Eg(„ (e„,_i_t)], G* is equal to the smallest G > such 
that, for some < a; < 1, (x, G) is a solution of (fTsT l and 
(fT4l i. where now f{x, G) is replaced by f{x, G). 

Using a proof technique analogous to that of Proposition |4] 
it is easy to show that the stability bound is still given by ( fTSl l, 
where now A2 = Y^lZT '9n.^2 



(ris.fe) 



and 



A. 




_'"'"^y^° 2} ^ ^^^^ _ 2^ 

u=l 



u ■ 



J{k,ns,k) 



is the expected number of weight-2 codewords of an {ns,k) 
linear block code picked uniformly at random in the ensemble 
of all [us , k) linear block codes without idle bits and with 
minimum distance at least 2. 

V. Numerical Threshold Optimization and 
Comparison with IRSA 

The analysis tool developed in Section|III]allows to calculate 
the threshold for a given choice of the Uc linear block codes 
"loh, h G {1, . . . , ric}, and of the p.m.f. P. Analogously, the 
tool developed in Section |IV] allows to evaluate the threshold 
of a CSA scheme under the random code hypothesis, for a 
given choice of the Smax lengths ris, n e {1, . . . , Smax}, and 
of the p.m.f. Q. These tools can be exploited to derive optimal 
(in the sense of maximizing the threshold G*) probability 
distributions P and Q in the two cases. 

Some optimized distribution profiles are shown in Table U 
Among the several possible algorithms available to find the 
global maximum of a nonlinear function, differential evolution 
[[191 has been used. In the upper part of the table, optimized 
probability distributions P are reported for an IRSA scheme 
with rates 1/2, 2/5 and 1/3, while in the lower part optimized 
probabilities distributions Q are illustrated for a CSA scheme 
with k = 2 and with the same rates, with the inclusion of 
R = 3/5, under the random code hypothesis. All distributions 
have been optimized under the constraint that the smallest 
local rate allowed for each user is 1/6. For each distribution, 
the corresponding threshold G* and stability bound (right- 
hand side of ( fTSl l) are shown. Note that in the CSA case the 
threshold values are average values: Specific choices of the 
codes 'i^h may lead to thresholds G* larger than G* . 
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Fig. 2. Thi'oughput versus the normalized offered traffic G for IRSA and 
CSA schemes with p.m.f.'s in Table |I] The bursts of each CSA user are split 
into fc = 2 segments, so that A^cSA = 2NsA- A^CSA = 1000, A^sA = 500. 



From Table H] we see that CSA is capable to achieve better 
thresholds than IRSA for i? = 1/2 and i? = 2/5, while for 
the lowest rate R = 1/3 IRSA exhibits a better threshold. 
Accordingly, the IRSA scheme seems to be preferable in the 
case of low rates R (i.e., for higher values of the excess 
power AP), while CSA is more interesting for higher values 
of R (i.e., when a higher power efficiency is required). Note 
also that CSA allows to achieve values of the overall rate 
R > 1/2, whereas only low rates R < 1/2 can be obtained 
from IRSA, unless some users transmit their burst in the MAC 
frame with no repetition. (In this latter case, however, no 
successful iterative IC can be guaranteed, so that we always 
have G* = 0.) For example, in Table J] an optimized CSA 
distribution of rate i? = 3/5 is reported. This distribution has 
no IRSA counterpart. 

To validate our asymptotic analysis, we performed nu- 
merical simulations in the case of a finite number M of 
users@ In Fig. [21 the throughput curves of IRSA and CSA 
schemes with the probability profiles from Table J] are de- 
picted as functions of the normalized offered traffic G. The 
throughput achieved by SA, S = Ge^'^ , is also shown 
for reference. Note that in our simulations for the CSA 
case, we combined the p.m.f.'s Q derived under the ran- 
dom code hypothesis with a specific choice of the com- 
ponent codes. In particular, we used linear block codes 
with the following generator matrices; G(3.2) = [110,011], 



G(4,2) = [1100,0111], G(5.2) = [11100,00111], G(8.2) = 

(9,2) = [111110000,011111111], 



[111100000,0111111], G 



G(i2,2) = [111111110000,0000011111111]. To stay fair, we 
compared CSA {k = 2) and IRSA schemes for the same 
frame duration Tp, which implies a number of slots A^csA 
twice the number of slots iVgA- Specifically, the simulations 
are for Nqsa = 1000 and A^sA = 500. For each value of G, 

^Here, one shall consider that each segment has to be encoded via a physical 
layer error correcting code before transmission on the MAC channel, and that 
the physical layer code for CSA is k times shorter than the corresponding 
code for IRSA. Thus, CSA may require working at slightly higher SNRs than 
IRSA, especially when short segments (and then short physical layer codes) 
are used. This aspect is not considered in this work. 



AI can be obtained from ([U. We observe a very good match 
between the asymptotic analysis and the simulations, the larger 
peak throughput of CSA than IRSA also for R — 1/3 being 
essentially due to the specific choice of the component codes 
(recall that G* is an average value). 

VI. Conclusion 

Coded slotted ALOHA has been introduced as a new 
opportunity for high-throughput random access to the MAC 
channel. Density evolution equations for CSA have been 
derived, optimal CSA schemes designed for several rates and 
their performance for a finite number of users simulated. The 
new scheme is promising when power efficiency is required. 
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